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00 ' Abstract 

The synchrony and variabihty have been discussed of the coupled Langevin 
^' model subjected to spatially correlated additive and multiplicative noise. We 

PJ , have employed numerical simulations and the analytical augmented-moment 

c/3 \ method which is the second-order moment method for local and global vari- 

ables [H. Hasegawa, Phys. Rev. E 67, 041903 (2003)]. It has been shown that 
the synchrony of an ensemble is increased (decreased) by a positive (negative) 
spatial correlation in both additive and multiplicative noise. Although the 
variability for local fluctuations is almost insensitive to spatial correlations, 
'^ ■ that for global fluctuations is increased (decreased) by positive (negative) cor- 

C \ relations. When a pulse input is applied, the synchrony is increased for the 

Q I correlated multiplicative noise, whereas it may be decreased for correlated 

additive noise coexisting with uncorrelated multiplicative noise. An applica- 
T^lj- . tion of our study to neuron ensembles has demonstrated the possibility that 

K^ \ information is conveyed by the variance and synchrony in input signals, which 

Lj. ' accounts for some neuronal experiments. 
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1 Introduction 

It has been realized that the coupled Langevin model is a valuable and useful model 
for a study of various phenomena observed in stochastic ensembles (for a recent 
review, see Ref. [1]). Independent (uncorrelated) additive and/or multiplicative 
noise has been widely adopted for theoretical analyses because of its mathematical 
simplicity. In natural phenomena, however, there exist some kinds of correlations in 
noise, such as spatial and temporal correlations, and a correlation between additive 
and multiplicative noise. In this paper we will pay our attention to the spatial 
correlation in noise. The Langevin model has been usually discussed with the use 
of the Fokker-Planck equation (FPE) for the probability distribution. In the case 
of correlated additive noise only, the probability distribution is expressed by the 
multivariate Gaussian probability with a covariance matrix. The effect of correlated 
additive noise has been extensively studied in neuroscience, where it is an important 
and essential problem to study the effect of correlations in noise and signals (for a 
review, see Ref. [2]). It has been shown that the synchrony and variability in 
neuron ensembles are much influenced by the spatial correlations [3]- [16]. The spatial 
correlation in additive noise enhances the synchrony of firings in a neuron ensemble, 
while it works to diminish beneficial roles of independent noise, such as the stochastic 
and coherent resonances and the population (pooling) effect [H UHl [HI [E] ; related 
discussions being given in Sec. 3. 

The problem becomes much difficult when multiplicative noise exists, for which 
the probability distribution generally becomes a non-Gaussian. Although an an- 
alytical expression of the stationary probability distribution for uncorrelated mul- 
tiplicative noise is available, that for correlated multiplicative noise has not been 
obtained yet. Indeed, only little theoretical study of the effect of spatially corre- 
lated multiplicative noise has been reported for subjects such as the noise-induced 
phase separation [17| and the Fisher information [l8]-[20], as far as the author is 
concerned. 

In a recent paper [21], we have studied stationary and dynamical properties of 
the coupled Langevin model subjected to uncorrelated additive and multiplicative 
noise. We employed the augmented moment method (AMM) which was developed 
for a study of stochastic systems with finite populations [22[ [23]. In the AMM, we 
consider global properties of ensembles, taking account of mean and fiuctuations 
(variances) of local and global variables. Although a calculation of the probability 
distribution for the spatially correlated multiplicative noise with the use of the FPE 
is very difficult as mentioned above, we may easily study its effects by using the 
AMM. It is the purpose of the present paper to apply the AMM to the coupled 
Langevin model including spatially correlated multiplicative noise and to study its 
effects on the synchrony and variability. 

The paper is organized as follows. In Sec. 2, we discuss the AMM for the spa- 
tially correlated Langevin model. With the use of analytical AMM and numerical 
methods, the synchrony and variability of the coupled Langevin model are inves- 
tigated. In Sec. 3, previous studies on the correlated multiplicative noise using 
the Gaussian approximation (18]-[20] are critically discussed. An application of our 



study to neuron ensembles is also presented with model calculations. The final Sec. 
4 is devoted to our conclusion. 



2 Formulation 

2.1 Adopted model 

We have assumed the A^-unit coupled Langevin model subjected to spatially corre- 
lated additive and multiplicative noise. The dynamics of a variable Xi {i = 1 to A^) 
is given by 

with 

«K^) = (I) E ^At)+m, (2) 

H{u) = -^==Q{u). (3) 

Here F{x) and G{x) are arbitrary functions of x, Z (= A^ — 1) denotes the coordi- 
nation number, /j(t) an input signal from external sources, w the coupling strength, 
and 0(m) the Heaviside function: 0(u) = 1 for u > and Q{u) = otherwise. 
We have included additive and multiplicative noise by C,i(t) and ^i(t), respectively, 
expressing zero-mean Gaussian white noise with correlations given by 

{v^{t)vAt')) = a%,+CM{l-S,,)]6{t-t'), (4) 

imUt')) = /3%,+CAil-6.,)Mt-t'), (5) 

{r]^{t) ^At')) = 0, (6) 

where the bracket (■) denotes the average, a (/3) expresses the magnitude of multi- 
phcative (additive) noise, and cm (ca) stands for the degree of the spatial correlation 
in multiplicative (additive) noise. Although our results to be present in the following 
are valid for any choice of H{x), we have adopted a simple analytic expression given 
by Eq. ([3]) in this study. 

We assume that external inputs have a variability defined by 

J,(t) = I{t) + 6h{t), (7) 

with 

(sm) = 0, (8) 

{6m6Ij{t')) = ^j[6,j + Si{l-6,j)]6{t-t'), (9) 

where 7/ and Sj denote the variance and degree of the spatial correlation, respec- 
tively, in external signals. We will investigate the response of the coupled Langevin 
model to correlated external inputs given by Eqs. ©-([H]). 



2.2 Augmented moment method 

In the AMM [221 [23], we define the three quantities of /i(t), 7(t), and p{t) expressed 

by 

p{t) = {x{t)) = j^Y.{^.{t)), (10) 

i 

lit) = ^T.il^^{t)-m?), (11) 

i 

Pit) = {[x{t)-p{t)r), (12) 

where X(t) = {l/N)J2iXi(t), l^if) expresses the mean, and 7(t) and pit) denote 
fluctuations in local (xj) and global variables (X), respectively. By using the Fokker- 
Planck equation (FPE), we obtain equations of motion for /i(t), 7(t) and p{t) which 
are given by (argument t is suppressed, details being given in the Appendix A) 

-^ = /o + /io + /27+ (yj [^0^1 + 3(^1^2 + ^0^3)7], (13) 

^ = 2h^ + ^iNp-^) + 2igl + 2gog2)a'^ + ^i + P' + a'gl (14) 



dp 

'di 



2/ip + 2hiwp + 2{gf + 2gog2)a'p 

1 Z 

+ — (7/ + /5' + a'(7o) + ^('5/7/ + CA/5' + CA/a%2), (15) 

where fi = {l/i\){d^F{p)/dx^), ge = il/i\)id^G{p)/dx^), h = {l/l\){d^H{u)/du^) 
and u = wp + /. Original A^-dimensional stochastic differential equations (DEs) 
given by Eqs. ([I])- ([3]) are transformed to the three-dimensional deterministic DEs 
given by Eqs. (TT3|) - (TT5|) . For •jj = Sj = ca = cm = 0, equations of motion given by 
Eqs. (IT51) - (TT^ reduce to those obtained in our previous study 
When we adopt F{x) and G{x) given by 





F{x) = -Xx, 


(16) 




G{x) = X, 


(17) 


Eqs. (T3D-(fT5) become 






dp 
dt 


9 
Xp + ho+ 2 , 


(18) 


d'y 


2hiwN ( 7 \ 9 


(19) 


dt 


^A7 -t- \p -t ^a ^ -\- Jr, 


dp 

dt 


(P + ZR) 
-2Xp + 2hiwp + 2a^ p+ ^ — ^, 


(20) 


with 


P = li + (^' + a'p\ 


(21) 




R = Sm + ca0^ + cmoi^P^, 


(22) 



where h^ = H{wp + I), and P and R express uncorrelated and correlated contribu- 
tions, respectively. We employ Eqs. (fT8l) - (!22l) in the remainder of this paper. 



2.3 Synchrony and variability 

2.3.1 Synchrony 

In order to quantitatively discuss the synchronization, we first consider the quantity 
S'(t) given by 

S'it) = 1^ E < [^^it) - ^At)? >= 2[7(t) - Pit)]. (23) 

When all neurons are in the same state: Xi{t) = X{t) for all i (the completely 
synchronous state), we obtain S'{t) = in Eq. fl23l) . On the contrary, in the 
asynchronous state where p = •y/N, it is given by S'{t) = 2(1 — 1/A^)7(t) = Sl^{t) 
We may define the normalized ratio for the synchrony given by 



Sit) ^l-^=(^)(P^-l-] (24) 

^'- S'.it) [z)[y{t) n)^ ^'^^ 

which is and 1 for completely asynchronous (5" = S'q) and synchronous states 
(S" = 0), respectively. 

2.3.2 Variability 

The local variability is conventionally given by 



'^''<'' - —^) — iw ^ ' 



where Sxiit) = Xj(t) — pit). Similarly, the global variability is defined by 



DAt) ^ ^^W^^^^CAt) 



Pit) 



(26) 



Pit) Pit) '''\At)' 

where 5X{t) = X{t) - pit). 

2.4 Stationary properties 

The stationary solution of Eqs. (fT8|l - (l20|) is given by 

^ = (A-tv2) ' ^'^^ 

^ P , ihw/Z)iP + ZR) 

^ 2(A - a2 + hiw/Z) 2(A - a^ + hiw/Z){X - a^ - hw) ' ^ ' 

_ P + ZR 

^ ~ 2N{X-a^-hw)' ^^^^ 

where pinP and R of Eqs. (EI]) and ([22]) is given by Eq. (E?]). We note in Eq. (E?]) 
that p is increased as I is increased with an enhancement factor of 1/(A — a^/2). A 



local fluctuation 7 is increased with increasing input fluctuations (7/) and/or noise 
(a, p) as Eq. ([28]) shows. In the limit oi Si = ca = cm = R = w = 0.0, Eqs. ([28]) 
and (129!) lead to p/7 = 1/A^ which expresses the central-limit theorem. From Eqs. 
and ([2SD, we obtain 



^ hwP + Z{\-a^)R 

P[Z{X^a^)-hiw{Z~l)] + hiwZR' ^ ' 



hiw 



Z{X-a^)-hiw{Z-l) 

Sill + caP"^ + CAf aV^ 
7/ + /52 + ^2^2 



for S"/ = Cyi = Cm = 0, (31) 

for w = 0, (32) 



where P and R in Eq. fl30]) are given by Eqs. fl2T]) and fl22]) . respectively. Equation 
( !30l) shows that the synchrony S is increased with increasing spatial correlations 
and/or the coupling. This is more clearly seen in the limit of no spatial correlations 
[Eq. ( 13T|) ] or no couplings [Eq. (l32l)]. The local and global variabilities, Cy and Dy, 
defined by Eqs. ( 125]) and ([26]) . respectively, are generally expressed in terms of P 
and R, and they are given for w = 0.0 by 

Cy = -I ,(;^_^.) j for^ = 0, (33) 



D 



y 




1 + ZSihi + (1 + Zca)/?' + (1 + ZcM)aV y^' 

iV(7/ + /5' + a^') / 

for w = 0. (34) 

The local variability Cy only weakly depends on the spatial correlation through the 
coupling, and it is independent of the correlation for w = 0. In contrast, the global 
variability Dy is increased (decreased) for positive (negative) correlations. In the 
limit oi Si = Ca = Cm = w = 0.0, Eq. fl51]) yields Dy = Cy/\/N expressing a 
smaller global variability in a larger- A^ ensemble (the population or pooling effect) 
[SliTOliniiTl]. 

The stability condition around the stationary state given by Eqs. (l27j) - (!29j) may 
be examined from eigenvalues of the Jacobian matrix of Eqs. (TT8]) - fl20]) . which are 
given by 

Ai = -X + — + hiw, (35) 

X2 = -2X + 2a'-^, (36) 

Zj 

A3 = -2A + 2a^ + 2/ii^. (37) 

The first eigenvalue of Ai arises from an equation of motion for /i, which is decoupled 
from the rest of variables. The stability condition for /x is given by 

hxw<{X-o?l2). (38) 



The stability condition for 7 and p is given by 

-Z{\-a^)<hiw <{\-a^). (39) 

Then for \ — a^ < hiw < A — a^/2, 7 and p are unstable but p remains stable. 
It is note that there is a limitation in a parameter value of c, as given by 

_ 1 < ( Sni + CA(3^ + CMa'^p'^ \ . . 

Z - [ 7, + /52 + a2^2 ) -^' l^uj 



which arises from the condition given by < p < 7 [see Eqs. ( 1T2l) and (l23i)]. 
When /3 = 7/ = 0, for example, a physically conceivable value of cm is given by 
-1/Z < Cm < 1. 

Cm dependence of the synchrony S is shown in Fig. [1] where a = 0.1, (3 = 0.1, 
c^ = 0.1, 7/ = S"/ = and A^ = 100. The condition given by Eq. (l40l) yields 
that —0.44 < Cm < 1.0 with /i = 0.5 for a given set of parameters. We note that 
the synchrony is increased with increasing s, and that the effect of the correlated 
variability is more considerable for larger p and w, as Eq. fl30|) shows. 

2.5 Dynamical properties 

In order to study the dynamical properties of our model given by Eqs. ([I])-®, we 
have performed direct simulations (DSs) by using the Heun method [211 l25j with a 
time step of 0.0001: DS results are averages of 100 trials. AMM calculations have 
been performed for Eqs. (TT8l) - (!20l) by using the fourth-order Runge-Kutta method 
with a time step of 0.01. We consider a set of typical parameters of A = 1.0, a = 0.1, 
P = 0.1, w = 0.5, 7/ = 5*7 = and A^ = 100. We apply a pulse input given by 

i{t) = Ae(t-4o)e(6o-t) + A, (41) 

with A = 0.4, Af, = 0.1, where 6(x) denotes the Heaviside function: B(x) = 1 for 
X > and otherwise. 

Figures [2](a)-l2](d) show time courses of p{t), •yit), S(t) and Cvit) for the corre- 
lated multiplicative noise (c^ = 0.0 and cm = 0.5). p{t) and S{t) are increased by 
an applied input at 40 < t < 60 shown by the chain curve in Fig. [2]^a), by which 
7(t) is slightly increased. The variability Cy(t) is decreased because of an increased 
p{t). The results of the AMM shown by the solid curves are in fairly good agreement 
with those of DS shown by the dashed curves. 

In contrast. Figs. [3](a)-l3](d) show time courses oi p(t), jit), S{t) and Cyit) for 
the correlated additive noise (ca = 0.1 and cm = 0.0). With an applied pulse input, 
pit) is increased and ■jit) is a little increased, as in the case of Figs. [2]|^a) and[2](b). 
However, the synchrony S{t) is decreased in Fig. [31(c) while it is increased in Fig. 
W^c). This difference arises from the fact that a decrease in S{t) in the former case 
is mainly due to an increase in P of the denominator of Eq. fl30l) , while in the latter 
case, its increase arises from an increase in R of the numerator of Eq. fl30l) . This 



point is more easily realized for w = 0, for which Eq. (!32l) yields 

S = -f^T^-T^, forc^ = 0, (42) 

""^^ for CM = 0. (43) 



/?2 + a2/i2 



The situation is almost the same even for finite w, as Figs. [2]^c) and [3](c) show. 
In both Figs. [2](d) and[3](d), Cyit) is decreased by an applied input because of an 
increased /i(t). 

3 Discussion 

3.1 A comparison with related studies 

We have investigated the stationary and dynamical properties of the spatially cor- 
related Langevin model given by Eqs. ([Il)-(l3]). In Ref. [26], we discussed the 
Fisher information in the Langevin model subjected to uncorrelated additive and 
multiplicative noise, which is a typical microscopic model showing the nonextensive 
behavior [27] . It is interesting to calculate the Fisher information of the Langevin 
model with correlated multiplicative noise. Such a calculation needs to solve the 
FPE of the Langevin model given by Eq. (lAip because the Fisher information is 
expressed in terms of derivatives of the probability distribution. For additive noise 
only (a = cm = 0), the stationary probability distribution p{{xk}) is expressed by 
the multivariate Gaussian distribution given by 



pi{xk}) oc exp 






(44) 



with the covariance matrix, Q, expressed by 

Qij = cr^[Sij + ca{1 - 6,j)], (45) 

where yU = H/X and a^ = j3^/2X. From Eqs. (jH]) and fH5|) . we obtain the Fisher 
information given by [18] 

9 = 2r. , ... — 7^^- (46) 

a^[l + [N -1}ca\ 

When multiplicative noise exists, a calculation of even stationary distribution 
becomes difficult, and it is generally not given by the Gaussian. The stationary 
distribution for uncorrelated additive and multiplicative noise {G{x) = x, ca = 
cm = 1i = 0.0) is given by [211 |26l [281 [29] 



In the limit of a = 0.0 and (3 ^ 0.0 {i.e. uncorrelated additive noise only), Eq. fj471) 
becomes the Gaussian distribution given by 

Pii^k}) oc l[e-^'/^"^^'^^-^^\ (48) 

i 

In the opposite limit of a 7^ 0.0, P = 0.0 and H > {i.e. uncorrelated multiplicative 
noise only), Eq. (1471) reduces to 

p{{x,}) oc n ^r^'^'^'^'^e-^''/"'" 0(^^), (49) 



yielding the Fisher information given by 



0"^ a^/i^ 



(50) 



where g = (2A + 3a^)/{2X + a^) and a| = aV/2A [26j. 

Unfortunately, we have not succeeded in obtaining the analytic expression for 
the stationary distribution of the Langevin model including corre/ateii multiplicative 
noise. In some previous studies [18]- [20], the stationary distribution for correlated 
multiplicative noise only {G{x) = x, ca = P = 0.0 and H = fi) is assumed to be 
expressed by the Gaussian distribution given by Eq. (jH]) with the covariance matrix 
given by 

Qij = (tIj iJiiJj[6ij + cm{1 - Sij)], (51) 

where /Xj (= (x,)) denotes the average of Xi and af^ a variance due to multiplicative 
noise. This is equivalent to assume that the multiplicative-noise term in the FPE 
given by Eq. (lAip is approximated as 

(y d d 

Y E T.[^iJ + Cm(1 - '^*j)]^^i^^J Pii^k}), 
Q/^ Q Q 

- Yimi'^»i + CM(l-'5*i)]^(a;»)^(xj)p({xfc}), 

Q/^ Q Q 

= yEE/^^/^it'^'J +^m(1 -'^u)] g^ g^ P({a^fc})- (52) 



If we adopt the Gaussian approximation given by Eq. fl52|) . with which multiplica- 
tive noise may be treated in the same way as additive noise, we obtain the AMM 
equations given by Eqs. ( IT8l) -( l20l) but without the third term of a^/i/2 in Eq. ( TTHj) . 
By using the Gaussian approximation given by Eq. ( ISTl) . Abbott and Dayan 
(AD) [T8] obtained the Fisher information expressed by [Eq. (4.7) of Ref. |18j ] 



NK 

alAl + {N - 1)cm] 

I 2N 

= -^T—^ — + ^ for A^^ 00, 54 

^ 2A^ r n 

= -^ — ^ H 5- for Cm = 0, (55) 



with K = N~^J2i[d'^nH{fii)/dfii]'^ = l//i^. Equation (!55|) is not in agreement with 
the exact expression given by Eq. (!50l) for uncorrelated multiphcative noise only. 

In stead of using the Langevin model, we may alternatively calculate the Fisher 
information of a spatially correlated nonextensive system by using the maximum- 
entropy method. In our recent paper [30] , we have obtained the analytic, stationary 
probability distribution which maximizes the Tsallis entropy |27j under the con- 
straints for a given set of the variance (a^) and covariance (c a'^). The Fisher 
information is expressed by [30] 

^ (56) 



a^[l + {N-l)cy 

= — - for AT ^ oo, (57) 

ca'' 

N 
= — fore = 0.0. (58) 

The Fisher information given by Eq. (|56ll is increased (decreased) by a negative 
(positive) correlation. This implies from the Cramer-Rao theorem that an unbiased 
estimate of fluctuations is improved by a negative spatial correlation, by which the 
synchrony is decreased as shown by Eqs. fl30|) and fl32l) . N and c dependences of 
the Fisher information given by Eq. fl56l) are different from those of qad given by 
Eq. fl53l) . although they are the same as those for additive noise only [Eq. fj46|) ]. 
It is noted that the Gaussian approximation given by Eq. fl^T]) or flS^ assumes 
the Gaussian distribution, although multiplicative noise generally yields the non- 
Gaussian distribution as shown by Eqs. (H7|) and (HOll . The spurious second term 
{2NK) in Eq. (1531) which is independent of cm and o"|,^, arises from an inappropriate 
Gaussian approximation. In discussing the Fisher information of spatially correlated 
nonextensive systems, we must take into account the detailed structure of the non- 
Gaussian distribution. 

3.2 Application to neuronal ensembles 

When 7/ and Sj in Eq. ^ are allowed to be time dependent, they may carry input 
information. This is easily realized if the AMM equations given by Eqs. flT51) - fl2Ul) 
are explicitly expressed in terms of /x, 7, and S as 

-f = -A/i + /io + — ^, 59 

at 2 

^ = _2A7 + 2/iiu;75 + 2a27 + 7j(t) + aV + /5', (60) 

at 

— = [7, (t) + « V + /5'] + - [7/ (t) Si (t) + CMa^iJ^ + CAf3^] 

at 7 7 

+ (H^)(l + ZS)(l-S). (61) 



which are derived with the use of Eq. \^£ 

10 



In order to numerically examine the possibility that input information is conveyed 
by ^i{t) and Si(t), we first apply a fluctuation-driven input given by 

7/(t) = fi0(t- 40)9(60 -t) + Eft, (62) 

with B = 0.4, Bb = 0.1, I{t) = 0.1 and Si{t) = 0.1 for A = 1.0, ca = 0.0, cm = 0.5, 
a = 0.1, P = 0.1 and A^ = 100. Time courses of /i(t), 7(t), S(t) and Cy(t) are shown 
in Figs. Il](a)-|l](d): the chain curves in Figs. ID^a) andlD^b) express I(t) and 7/(t), 
respectively. When the magnitude of 7/(t) is increased at 40 < t < 60, 7(t) and 
Cv{t) are much increased, while there is no changes in fi{t) because it is decoupled 
from the rest of variables in Eq. flTSj) . S{t) is slightly modified only at t ~ 40 and 
t ~ 60 where the 7/(t) is on and off. 

Next we apply a synchrony-driven input Si{t) given by 

Sj{t) =Ce{t- 40)9(60 -t) + Cb, (63) 

with C = 0.4, Cb = 0.1, I{I) = 0.1 and 7/(t) = 0.1. Figures [5|a)-Hd) show time 
courses of /i(t), ■j{t), S{t) and Cv(t): the chain curves in Figs. M^a) and[5]^c) express 
I{t) and Si{t), respectively. An increase in synchrony-driven input at 40 < t < 60 
induces a significant increase in S{t) and slight increases in •yit) and Cyit), but no 
changes in fi(t). 

When we regard a variable Xj in the Langevin model given by Eqs. ([I])- ([3]) as 
the firing rate r^ (> 0) of a neuron z in a neuron ensemble, our model expresses 
the neuronal model proposed in Refs. [211 [32] • It belongs to the firing-rate (rate- 
code) models such as the Wilson-Cowan [33] and Hopfield models [31], in which 
a neuron is regarded as a transducer from input rate signals to output rate ones. 
Alternative neuronal models are spiking-neuron (temporal-code) models such as the 
Hodgkin-Huxley [35] , FitzHugh-Nagumo [36l [37] and integrate-and-fire (IF) models 
[38] . Various attempts have been proposed to obtain the firing-rate model, starting 
from spiking-neuron models [39]- [l3]. It is difficult to analytically calculate the firing 
rate based on spiking-neuron models, except for the IF- type model [.38] . It has been 
shown with the use of the IF model that information transmission is possible by 
noise-coded signals [lH US], and that the modulation of the synchrony is possible 
without a change in firing rate [46] . Model calculations shown in Figs. H] and [S] 
have demonstrated the possibility that information may be conveyed by 7/(t) and 
Si{t), which is partly supported by results of the IF model [11]-[16]. Some relevant 
results have been reported in neuronal experiments [lT]-[52]. In motor tasks of 
monkey, firing rate and synchrony are considered to encode behavioral events and 
cognitive events, respectively [1?]. During visual tasks, rate and synchrony are 
suggested to encode task-related signals and expectation, respectively [IB]- A change 
in synchrony may amplify behaviorally relevant signals in V4 of monkey [l9]. The 
synchrony is modified without a change in firing rate in some experiments [171 UHl 
[51] . Synchrony-dependent firing-rate signal is shown to propagate in iteratively 
constructed networks in vitro 
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4 Conclusion 

With the use of DSs and the AMM [22l [23] , the effects of spatially correlated addi- 
tive and multiplicative noise have been discussed on the synchrony and variability 
in the coupled Langevin model. Our calculations have shown the following: (i) the 
synchrony is increased (decreased) by the positive (negative) correlation in additive 
and multiplicative noise [Eq. (!30|) ]. (ii) although an applied pulse input works to in- 
crease the synchrony for correlated multiplicative noise, it is possible to decrease the 
synchrony when correlated additive noise coexists with uncorrelated multiplicative 
one, (iii) the local variability Cy is almost independent of spatial correlations, while 
global variability Dy is increased (decreased) with increasing the positive (negative) 
correlation, and (iv) information may be carried by variance and synchrony in input 
signals. The item (iv) is consistent with the results of Refs. [11]-[16] and elucidates 
some phenomena observed in neuronal experiments |47j-|52j. 

Although we have applied the AMM to the Langevin model in this paper, it 
is possible to apply it to other types of stochastic neuronal models such as the 
FitzHugh and Hodgkin-Huxley models subjected to correlated additive and multi- 
plicative noise, which is left as our future study. 
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Appendix A: Derivation of the AMM equations 

The Fokker-Planck equation (FPE) for the Langevin equation given by Eqs. 
dl])-© in the Stratonovich representation is expressed by pTf [53] 

-^^p{{xk],t) = -Y,—{[F{x,) + H{u.;)]p{{xk},t)} 
+ jT.J: d^M^^^ + ^^(1 - ^«^)] P(M^ t)} 

1 J J 

+YEE^{[^^. + CA(i-5.,)]p(W,t)} 

i j 'I- i 

where Ui = {w/Z) J2j(j^i) ^j + ^ and H'{ui) is absorbed in a new definition of 77- in 
its second term. 

Equations of motion for moments, (xj) and (xj Xj), are derived with the use of 
FPE [21]: 

-^ = {F{x.) + H{u.)) + —{G\x.)G{x.)), (A2) 
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+ -^\kx^G\x,)G{x,)) + i,x,G\x,)G{x,))\ 

+ (l-(5,,)['^/7/ + CA/5' + CAf«'(G(x,)G(x,))]. (A3) 

In the AMM [221 123], the three quantities of /x, 7 and p are defined by Eqs. 
(ITU1) - (TT^ . We use the expansion given by 

Xi = n + Sxi, (A4) 

and the relations given by 

d/i _ 1 ^d{ri) 

dl _ 1 v^ d{{5rif) 

dp _ 1 (U^TjSrj) 

^ - Jf2l.l^^^^- (A7) 

For example, Eq. (jASP for dp/dt is calculated as follows: 

l^(F(r,)) = /0 + /27, (A8) 

]^Y.^H{u,)) = ho, (A9) 

^$:(G"(r,)G(r,)) = ^0^1 + 3(^70^3 + ^1^2)7- (AlO) 

Equations ( 1A6P and ( 1A7P are calculated in a similar way. Then, we have obtained 
equations of motion for p(t), ■j(t) and p(t) given by Eqs. (TT3l)-(TT5l). 
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Figure 1: (Color online) Cm dependence of the stationary synchrony S ioi w = 0.5 
(solid curves) and w = 0.0 (dashed curves) with a = 0.1, /3 = 0.1, ca = 0.1, 
ji = Si = and N = 100, fj, being treated as a parameter. 



Figure 2: (Color online) Time courses of (a) fi{t), (b) 7(t), (c) S{t) and (d) Cv(t) 
with the correlated multiplicative noise {ca = 0.0, cm = 0.5, a = 0.1, /3 = 0.1) for a 
pulse input given by Eq. (HTj) with A = 0.4, Af, = 0.1: the solid and dotted curves 
express results of the AMM and DS, respectively: the chain curve in (a) expresses 
an input of /(t) (A = 1.0, ^^ = 7^ = 0.0 and N = 100). 



Figure 3: (Color online) Time courses of (a) /i(t), (b) 7(t), (c) S(t) and (d) Cyit) 
with the correlated additive noise (ca = 0.1, cm = 0.0, a = 0.1, P = 0.1) for a pulse 
input given by Eq. (HT!) with A = 0.4 and A^ = 0.1: the sohd and dashed curves 
denote results of the AMM and DSs, respectively: the chain curve in (a) expresses 
an input of I{t) (A = 1.0, ^^ = 7^ = 0.0 and A^ = 100). 



Figure 4: (Color online) Time courses of (a) /i(t), (b) ■y(t), (c) S(t) and (d) Cy(t) 
for a fluctuation-driven input of ■Jiit) given by Eq. (!62l) with Sj = 0.1, B = 
0.4 and Bf, = 0.1: the solid and dotted curves express results of the AMM and 
DS, respectively: the chain curves in (a) and (b) express inputs of I(t) and 7/(t), 
respectively (A = 1.0, ca = 0.0, cm = 0.5, a = 0.1, (3 = 0.1 and A^ = 100). 



Figure 5: (Color online) Time courses of (a) /i(t), (b) 7(t), (c) S(t) and (d) Cyit) for 
a synchrony-driven input of Sj(t) given by Eq. (1^ with 7/ = 0.1, C = 0.4 and Cb = 
0.1: the solid and dotted curves express results of the AMM and DS, respectively: 
the chain curves in (a) and (c) express inputs of /(t) and Si{t), respectively (A = 1.0, 
CA = 0.0, Cm = 0.5, a = 0.1, p = 0.1 and A^ = 100). 
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